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Exact scheme independence at two loops 
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We further develop an algorithmic and diagrammatic computational framework for very general 
exact renormalization groups, where the embedded regularisation scheme, parametrised by a general 
cutoff function and infinitely many higher point vertices, is left unspecified. Calculations proceed 
iteratively, by integrating by parts with respect to the effective cutoff, thus introducing effective 
propagators, and differentials of vertices that can be expanded using the flow equations; many 
cancellations occur on using the fact that the effective propagator is the inverse of the classical 
Wilsonian two-point vertex. We demonstrate the power of these methods by computing the beta 
function up to two loops in massless four dimensional scalar field theory, obtaining the expected 
universal coefficients, independent of the details of the regularisation scheme. 
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I. INTRODUCTION AND CONCLUSIONS 

The deeper understanding of renormalization, due to 
Wilson, follows most directly in the continuum from 
the exact renormalization group (ERG) flow equations 
0. The fact that solutions of these equations, for 
the Wilsonian effective action S, can be found directly 
in terms of renormalised quantities, that all physics 
(e.g. Green functions) can be extracted from S, and 
that renormalizability is trivially preserved in almost 
any approximation 0, Q , turns these ideas into a pow- 
erful framework for considering both perturbative and 
non-perturbative approximations (see for ex amp le refs. 

In the past a number of different versions and ways of 
deriving the ERG have been proposed [JEHSEIlHH 
IT(i l2l| , which however have been shown to be equivalent 
under changes of variables |2|, y> LLll LUJ LUj ■ Recently, far 
more general versions of the ERG have been considered 
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(13| . All the ERGs, including these generalised ones, can 
be seen to be parametrised by a functional \& , that 

induces a rcparamctrisation (field redefinition) along the 
flow, and acts as a connection between the theory space 
of actions at different effective cutoff scales A. As a re- 
sult, local to some generic point A on the flow, all these 
ERGs may be shown to be just reparametrisations of 
each other. When this reparametrisation can be extended 
globally, the result is an immediate proof of scheme inde- 
pendence for physical observables. Indeed computations 
of physical quantities then differ only through some field 
reparametrisation. 

One practical example is an explicit field redefinition 
that interpolates between results computed using differ- 
ent choices of cutoff function [l^. Even more dramatic 
than this however, is the use of this freedom to adapt 
the ERG to certain forms of approximation or special 
physical problems 0. In particular, recently there has 
been substantial progress in adapting these ideas to gauge 
theory. It turns out that not only can one introduce 
an effective cutoff A in a way that does not break the 
gauge invariance but careful choices of 4" allow the 
gauge invariance to be preserved manifestly, in fact not 
even gauge fixed, along the flow and in the solutions S 

mum 

Ref. 01 did not answer the question of precisely when 
the automatic local equivalence of two ERGs can be ex- 
tended globally. In ref. ^(|, we showed that for four 
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dimensional one-component scalar field theory, the uni- 
versal one-loop (3 function is obtained for a general form 
of 'J, involving a very general 'seed' action S. For sim- 
plicity we chose to keep ip <-» — ip invariance (p being the 
scalar field) and the S two-point vertex was specified to 
be equal to that of the classical effective action, which 
thus determines them, up to a choice of cutoff function 
c. The only further requirements we imposed were that 
the vertices of S be infinitely differentiable and lead to 
convergent momentum integrals. These easily met re- 
quirements are needed in any case for the ERG equation 
to make sense at the quantum level [lq|. 

In this paper we will show that the universal two-loop (3 
function also comes out correct for this very general form 
of \&. It is natural to conjecture then, that to all orders 
in perturbation theory, the only constraints required on 
S to get universal answers for physical quantities are the 
ones in italics above. Furthermore, it is surely possible 
to show in perturbation theory, that all such ERGs are 
reparametrisations of each other, and indeed to construct 
the map perturbatively. 

We will see that the iterative diagrammatic compu- 
tational framework introduced in ref. ^(|, and further 
refined in the gauge theory context in ref. |15| . extends 
straightforwardly to higher loops. The method works by 
turning the large redundancy in (here encapsulated in 
S) to our advantage. Since we are not allowed to inquire 
into the form of its vertices, the calculational steps are 
severely limited. The inherent generality thus acts as a 
roadmap to the most efficient computation. Indeed, since 
the form of the vertices is never specified, the majority 
of the calculation is best performed by manipulating the 
diagrams themselves. 

However, it should be emphasised that, although the 
calculational procedure is the most efficient, the actual 
computation of the two-loop /3-function, presented here, 
could be considerably shorter. Indeed, if our purpose 
were just to compute the two-loop /^-function then, hav- 
ing used the redundancy of if? to uncover the best calcu- 
lational procedure, we could use the simplest form of S 
sufficient to yield a valid ERG. This corresponds to dis- 
carding all seed action vertices, other than the classical 
two-point vertex, reducing the ERG equation to Polchin- 
ski's version. We would then be left with a powerful al- 
gorithmic and diagrammatic procedure, which does not 
require the cutoff function c to be specified. 

Our aim is rather to explicitly demonstrate scheme in- 
dependence and so we do not restrict S in this way. In- 
deed, we even allow S to have its own completely unspec- 
ified (loop) expansion in h. Furthermore, by retaining a 
generic seed action, we gain valuable experience for the 
gauge theory case, where S necessarily contains many 
interaction vertices [l5| . 

An important prerequisite for this method to work, is 
that the effective action is written in self-similar form 
p??! , in which no other explicit scale is introduced except 
A. (Note that although we formulate the method for a 
massless theory, nothing precludes the method from be- 



ing applied to massive theories. The required self similar 
flow is achieved in this case simply by introducing the ap- 
propriate RG concept of a running mass, i.e. a function 
of A with its own (3 function, see for example @,EJ-) 

We also constrain the form of the flow equation so that 
in this context, the kernel appearing in the ERG equa- 
tion, is the differential of the effective propagator, the 
latter being the inverse of the classical effective action 
two-point vertex. (In gauge theory this condition holds 
only up to gauge transformations 0,0.) 

The computation then proceeds as follows. We intro- 
duce effective propagators in diagrams containing no ex- 
plicit seed action vertices, by integrating by parts with 
respect to A. This results in total A-derivative contri- 
butions, and terms where the A-derivative acts on effec- 
tive action vertices. We use the flow equation to process 
these latter further. Then, using the fact that the effec- 
tive propagator is the inverse of the classical two-point 
vertex, many cancellations typically occur. The above 
procedure is repeated until there is no explicit depen- 
dence left upon the (generic) seed action S. The whole 
of this iterative procedure may be performed entirely di- 
agrammatically. 

Universal terms arise from the total A-derivative con- 
tributions. Whenever the total A-derivative acts on con- 
vergent momentum integrals, the result follows trivially 
by dimensions; in particular, we frequently exploit the 
fact that dimensionless cases simply vanish. Although ul- 
traviolet finiteness is built in to the ERG, non- vanishing 
universal terms arise from dimensionless momentum in- 
tegrals that are infrared convergent only after the A- 
derivative is taken. At the two-loop level, we generate for 
the first time, contributions that are infrared divergent 
even after differentiation with respect to A. Of course at 
finite A, everything is infrared finite in the ERG, again 
by construction. Thus the game here is to rearrange so 
these infrared divergences cancel amongst each other; in- 
telligent combinations of these terms then are convergent 
and vanish, or at worst result in calculable universal con- 
tributions. 

In this way, the standard two-loop (3 function coeffi- 
cient, — 17/3(47r) 4 , is derived, however without specifying 
at any stage the form of the cutoff function or the form 
of any of the higher point vertices. 

We note that this two-loop (3 function coefficient has 
already been derived within the ERG, with various cutoff 
functions, and in various ways, corresponding to differ- 
ing motivations [H E3, El E3 HI Ej Although ref. 
|2ll | considered a general cutoff function, this work is the 
first to treat the more general case of an arbitrary seed 
action S (corresponding to a continuum version of an ar- 
bitrary "blocking scheme" 0), and to reduce the compu- 
tation in this general framework to a largely algorithmic 
and diagrammatic approach. We expect these insights to 
be especially useful for higher loop computations in the 
gauge invariant ERGs [25j. 



3 



II. A GENERALISED EXACT 
RENORMALIZATION GROUP EQUATION 

We will consider a massless one-component scalar field 
theory in four Euclidean dimensions, and sketch the 
derivation of a generalised ERG equation, starting from 
Polchinski's equation Q. For details we refer the reader 
to 0. 

The Wilsonian renormalization group (RG) is defined 
in terms of some effective ultraviolet cutoff A [ij . Polchin- 
ski's version of Wilson's ERG equation [l| implements 
this transparently through a cutoff function c(p 2 /A 2 ), 
which modifies propagators 1/p 2 to A = c/p 2 . We take 
c to be smooth, i.e. infinitely differentiable, and always 
positive. It satisfies c(0) = 1 so that low energies are un- 
altered, and tends to zero as p 2 /A 2 ~ * oo sufficiently fast 
that all Feynman diagrams are ultra-violet regulated. 

The partition function is given as the functional in- 
tegral of the measure, e~ s , where S is the Wilsonian 
effective action 



%;A] = i/ 



d 4 P 



pVV + ^foA], (1) 



Cp 1 = c (p 2 /A 2 ). The first term in the above, namely 
the regularised kinetic term, will be referred to as the 
seed action and denoted by S. 

Demanding that physics be invariant under the renor- 
malization group transformation Ah> A - SA, results in 
a functional differential equation for the effective inter- 
action , which can be recast in terms of the total effec- 
tive action, S = S + S mt [cf. eq. ©], E = S - 2S, and 
the differentiated effective propagator A = — A<9aA = 
£c>(p 2 /A 2 ), as 

S = -Ad S--— A — -!— A — (2) 

2 Sip Sip 2 Sip Sip 

(up to a vacuum energy term discarded in 0. We will 
also refer to A as the kernel.) 

In the above, prime denotes differentiation with respect 
to the function's argument (herep 2 /A 2 ) and the following 
shorthand has been introduced: for any two functions 
f(x) and g(y) and a momentum space kernel W(p 2 /A 2 ), 
with A being the effective cutoff, 



f-W-g 



a^xdSj f(x)W xy g(y), 



(3) 



where W xy = f-£^W(p 2 /A 2 )e ip < x - y \ (S may therefore 

be written as h d^ip-c^ 1 -d^tp.) 

Eq. @, which is given diagrammatically as in fig. Q 
can be turned into a flow equation for the measure, 



An s IS if <5E a 
A9A6 =-2^' A 'U e 



(4) 





which makes the invariance of the partition function un- 
der the renormalization group transformation manifest 



FIG. 1: Graphical representation of the ERG equation. S, E 
ra-point vertices are represented by circles labelled respectively 
by S, £ with n legs attached, while bulleted lines stand for A. 



by showing that the measure flows into a total functional 
derivative. 

Eqs. J2J, |@J may also be reinterpreted in a slightly 
different way: introducing W — — A |^ , we can 
regard the change in the partition function as correspond- 
ing to the field redefinition ip — > ip + 5A ty, from which we 
can infer that integrating out degrees of freedom is just 
equivalent to a reparametrisation of the partition func- 
tion |l3t IT^j . This somewhat counterintuitive result is 
allowed in the continuum because of the infinite number 
of degrees of freedom per unit volume. 

Different forms of ERG equations are obtained by 
choosing different '5, and there is a great deal of freedom 
in such a choice of "blocking scheme" 01 . Nonetheless, 
physical quantities should not depend on the particular 
scheme, in other words they should turn out to be uni- 
versal. Therefore we can choose a <J/ which is best suited 
for our purposes (for example one that generates a mani- 
festly gauge invariant ERG ^3,^3) ana ^ s ^ t ne right 
answer when computing physical quantities. 

Moreover, if physical results are to come out the same, 
irrespective of the choice of 4", we might well decide not to 
pick a specific form for it, but rather leave it general (ex- 
cept for some very general requirements discussed later) 
and implement a method for calculating universal quan- 
tities which does not depend on the details put in by 
hand. The main advantage of such a procedure is that 
we must be able to see that all the dependence of our 
results on unspecified quantities clearly cancels out. 

In this paper we want to go along these lines and gener- 
alise our previous results 16] by calculating the two-loop 
beta function in a scalar field theory with a generalised 
'J. With so much freedom we have to restrict it to be 
able to be concrete; we choose to consider a general seed 
action S of the form outlined below. 

Firstly, we want the tree-level two-point vertex in S 
to be the same as the required two-point vertex for the 
classical effective action. Since we are dealing with a 
massless theory, that means that both two-point vertices 
should be set to be equal to the regularised kinetic term 
in eq. 0^2.9). We will see in the next section that this 
constraint can be consistently imposed on solutions of 
the flow equation J5J. 

(In general we impose that the two-point vertices of 
S and the effective action must coincide at the classical 
level. The form of the flow equations then imply that A 
really is the differentiated propagator in this context, i.e. 
the differential of A, where A is the inverse of the clas- 
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sical effective action two-point vertex. Here we aiready 
soived this constraint in writing eq. @- However, for a 
massive theory, any classical mass would have to appear 
already in A, thus slightly generalising eq. Further- 
more, in manifestly gauge invariant ERGs |l4l Il5j . the 
inverse of the classical two-point vertex does not exist 
for the gauge fields; instead, A is the inverse, only up 
to a gauge transformation. It is not necessary that the 
two-point vertices of S and the classical effective action 
coincide. We choose to require this purely for the sig- 
nificant technical advantages it brings to our method, as 
described later.) 

Secondly, we choose to leave the tp *-* — <p symmetry 
alone, which implies that S must be even under this sym- 
metry. We are left with a generalised exact renormaliza- 
tion group parametrised by the infinite set of seed action 
2n-point vertices, n > 2. We will leave each of these 
vertices as completely unspecified functions of their mo- 
menta except for the very general requirements that the 
vertices be infinitely differentiable and lead to convergent 
momentum integrals. The first condition ensures that no 
spurious infrared singularities are introduced and that 
all effective vertices can be Taylor expanded in their mo- 
menta to any order 0,0]. The second condition is neces- 
sary for the flow equation to make sense at the quantum 
level and also ensures that the flow actually corresponds 
to integrating out modes |13l Il4j . 

We are therefore incorporating in the momentum de- 
pendence of each of the seed action 2n-point vertices an 
infinite number of parameters. Of course these infinite 
number of vertices, each with an infinite number of pa- 
rameters, then appear in the effective action S as a con- 
sequence of the flow equation. Remarkably, we can still 
compute the two-loop (3 function. 

III. PRELIMINARIES 

A. Self-similar flow 

We denote the vertices of the effective action as 

S^(p;A) = S^( Pl ,p 2 ,--- ,p 2n ;A), (5) 

where we have factored out the momentum conserving 5 
function: 

(6) 

(and similarly for S). 

In a standard perturbative treatment, we would de- 
fine the theory, by stating that at the classical level the 
Lagrangian density takes the form 

\{d,p) 2 + ^. (7) 

After regularising the theory, the coupling is replaced by 
a bare coupling, and a bare mass is introduced, and these 



are solved for, order by order in A, so that at physical 
scales all quantities of physical interest {e.g. Green func- 
tions) are finite and correspond to a massless theory. An 
essential step in this treatment is to define what we really 
mean by A through some renormalization condition. 

If we regularise the theory, by using the cutoff function 
c at some ultra-violet cutoff scale A = Ao, then the Wilso- 
nian effective action at physical scales S[p; A] will simply 
result from integrating eq. J2J), with the initial condition 
that S[(p; Ao] is the bare action. This can be done order 
by order in A. In the continuum limit Ao — > oo, we obtain 
the Wilsonian effective action expressed in renormalized 
terms. 

We will not approach the calculation this way. Instead, 
firstly we choose renormalization conditions that only in- 
volve the scale A, for example in this case we use 

S^(p; A) ee S (2 Hp, -p; A) = ct(A)A 2 + p 2 + 0(//A 2 ), 

(8) 

S (4) (0;A) = A, (9) 

where a — o~\ A + a 2 A 2 + ■ ■ • is a function of A which we 
determine self-consistently. This sets the physical mass 
to zero implicitly by ensuring that the only scale that 
appears is A, which itself tends to zero as all momenta 
are integrated out [3(j • In order to obtain unit p 2 coeffi- 
cient in eq. (JSJ), of course we also have to reparametrise 
the field by the wavefunction renormalization factor. In 
the limit that Ao — > 00, apart from the A dependence ex- 
pected by naive (a.k.a. engineering) dimensions, the only 
dependence on A then appears through A (A). Equiva- 
lently, if we rescale the field and positions/momenta by 
the appropriate (engineering) powers of A to make ev- 
erything dimensionless, then the effective action takes 
the self-similar 27] form jS[^j;A], in which the only de- 
pendence on A is through the coupling A (A). The ex- 
istence of such a solution is equivalent to a statement 
of renormalizability, and the solution corresponds to the 
renormalized trajectory [3l| 

Secondly, we considerably simplify the analysis, both 
conceptually and in terms of procedure, by solving the 
ERG directly in terms of the self-similar solution S[tp; A], 

order by order inA@,H[ll0HH[ll[l3. In ^ 

way, we never introduce an overall cutoff (e.g. Aq) or the 

notion of a bare action at all. Instead, we obtain the 
continuum physics directly. 



B. Perturbative expansion 

Rather than directly use the form of the flow equation 
in eq. ©, we first make explicit the wavefunction renor- 
malization contribution. Moreover, we rescale ip so as to 
put the coupling constant in front of the action. This will 
ensure the expansion in the coupling constant coincides 
with the one in ft, the actual expansion parameter being 
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in fact \h. Our flow equation then reads 0] 



-Ad A 



A 



j_ 

2A 



1 tf(5-2S) ^ 55 _ 1_6_ • 6(S-2S) 



(10) 



■4 

where J p = J,^ ? 4 and a tilde has been put over rescaled 



as 



quantities — so for example ip = -j^p arid S 1 
well as over the anomalous dimension 7 (7 = A<9a log Z , 
Z being the wavefunction renormalization) to signify it 
has absorbed the change to A<9a 1 5 ■ 

(Note that eq. is actually not the result of rescal- 
ing the wavefunction renormalization out of eq. © . We 
must in addition change the cutoff function c 1— » cZ in 
the flow equation 0,0, E|. However, the important 
point is that it clearly still satisfies the requirements of 
an ERG equation: namely that the partition function is 
invariant under the flow and that the flow corresponds to 
integrating out degrees of freedom [T3L Il5| . Hence it is a 
valid and even more appropriate starting point when the 
wavefunction renormalization is to be taken into account. 
This is another example of the immense freedom in the 
choice of the ERG equation.) 

In order to simplify the notation, the tildes will be 
removed from now on. In these new, rescaled variables, 
the renormalization conditions are set to 

S {2) ( P - A) = S {2) ( P , -p; A) = cr(A)A 2 + p 2 + 0(p 4 /A 2 ), 

(11) 



S (4) (0;A) = 1. 



(12) 



Both conditions are already saturated at tree level. (To 
see this it is sufficient to note that, since the theory is 
massless, the only scale involved is A. Since the tree-level 
four-point vertex is dimensionless, it must be a constant 



at null momenta. Thus S^ 4 -* (0; A) 



4 4) (0;A ) 



1, 



where the lower index signifies the coupling is intended 
at the tree level. Similar arguments apply to the tree- 
level two-point function.) 

Expanding the action, the seed action, the beta func- 
tion /3(A) = Ac?aA and the anomalous dimension in pow- 
ers of the coupling constant, and bearing in mind that the 
(inverse of the) coupling constant now appears in front 
of S: 



S[p;X\ = S + XS 1 + \ 2 S 2 + 

S[ip;X] = S + XS 1 + \ 2 S 2 + 

/3(A) = /3 1 A 2 +/3 2 A 3 + -- - , 

7(A) = 7iA + 7 2 A 2 H 



yields the loopwise expansion of the flow equation 



02 + P2S0 



71 

2 (p ' Sip 



So = 



1 SSo 1 ST, 

2 Sip Sip ' 



Si + PiSq + y "f~^ 



72 SSo_ 

2 (p ' Sip 



SSi 
Sip 
5S 2 
Sip 



•A- 



•A- 



S_Uo 


_SSo 


A 


SSi 


1 s 


Sip 


Sip 


Sip 


2S^p" 


Slip 


^SSo 


A 


ss 2 


1SS! 


Sip 


Sip 


Sip 


+ 2~Slp 



SZp 

Sip ' 
(5Ei 
Sip 



1 5 



ST,! 



2 Sip Sip 



(13) 
(14) 
(15) 



etc. where again dots above quantities signify — A<9a, 
S rl = Sn 2Sn and n n = Sn Sn* 

(3 and 7, at one- and two-loop order, may be extracted 
directly from eqs. I|14fl . (|15fl . as specialised to the two- and 
four-point effective couplings, once the renormalization 
conditions have been taken into account. 

The procedure is very straightforward: as the renor- 
malization conditions are already saturated at tree level, 
there must be no quantum corrections to the two-point 
effective coupling at order p 2 , nor to the four-point at 



null momenta, i.e. 

S< 2 >(p;A) =S( 4 >(0;A) = Vn > 1, (16) 

where the notation f\ p2 signifies that the coefficient of 
p 2 in the series expansion of / must be taken. Hence 
the ERG equations for these special parts of the quan- 
tum corrections greatly simplify, reducing to algebraic 
equations which can be solved for /?, 7 at any order of 
perturbation theory. 

As explained in the previous section, our S is not com- 
pletely arbitrary. Apart from the very general require- 
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merits on the differentiability and integrability of its ver- 
tices mentioned earlier, for convenience we restrict S to 
have only even-point vertices, and constrain its tree-level 
two-point vertex so that it is equal to that of the classical 
effective action: 

S^(p) = S^(p). (17) 

Although this constraint is not necessary, it greatly sim- 
plifies the flow equations for higher point vertices, as it 

(2) 

implies, for example, that IIq = 0. Recall that we also 
wish to set 



(18) 



(2) 

However, Sq (p) is already determined, up to an integra- 
tion constant, by the flow equation. Using eq. (|17fl and 
the two-point part of eq. (|13f) . and rearranging, we have 



(19) 



which shows that eq. (Jffi is consistent with the flow equa- 
tion. Indeed, 

4 2) (p)A p = l, (20) 
a central relation in the calculation that follows. 



1 




So 



FIG. 2: Diagrammatic representation of eq. 1201 . 

Before going to the details of the two-loop calculation, 
to which the next section will be devoted, we will describe 
what our method consists of by rederiving the one-loop 
contribution to /3 0|. Being a much simpler calculation, 
it constitutes the perfect ground for illustrating our strat- 
egy and, moreover, it will help the reader to get familiar 
with the diagrammatics. 

In order to ease notation in the figures, in what follows 
the vertices of the effective action will be labelled by their 
loop order only, and those of S by their loop order with 
a "hat" on top. 



where A = ^c'(O). Eqs. J3T}, are represented 

diagrammatically as in figs. [210] 



0i + 27i 





FIG. 3: Graphical representation of eq. (I2H . All the external 
legs have null momentum. 



A + 71 = -2 -H i 



FIG. 4: Graphical representation of eq. 1221 . 

We start processing eq. Ij22(l by integrating by parts 
the diagram containing the four-point Sq vertex. 




ft. +7i 



1 

+ 2 



A„ 



xSo'{p,-p,q,-q) + A q S^ 4) [p, -p, q, -q) 



2S[ 2 \p)A p S^>(p) 



(2), 



(23) 






V 


_ 1 

2 




P2 





C. One-loop beta function 

Specialising eq. Ijl4|l to the two- and four-point effective 
couplings and imposing the renormalization conditions, 
eq. 116fl . yields 

ft + 2 7l = iS[ 2) (0) Ao 4 4) (0) - 4^ 4) (0) A S[ 2) (0) 



\J^A q ^\d, q ,- q ), 



(21) 



ft +7i 



2S^(p)A p Sr>( P ) 



(2), 



- \ J Ay Sq 4) (p, ~P,q,-q) 



(22) 



+ i 




P2 




p2 



FIG. 5: Integrating by parts the four-point vertex in eq. (1221) 
[cf. eq. {Sty. 

Now, the first term in the above represents a fully con- 
vergent integral [cf. eq. (I23[) ]. therefore the order of the 
derivative and integral signs can be exchanged. More- 
over, as the integrand is dimcnsionlcss, there can be no 
dependence upon A after the momentum integral is car- 
ried out. Hence the result vanishes identically. The sec- 
ond term can be processed further by making use of the 
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tree-level flow equation for Sq (p, — p, q, — g): 



- A q S ( ^\p, ~p,q,~q) 
?(2) 



x(-2S^(p)A p -2S^(q)A q 



(2)/ 



Scj 4) (P) -P) 9, -<?) 

= -A„ 



x -5'o 4) (P) _ P' <?' ~<?) - A 9 s£ 2) (p) A p S^ 4) (p, -p, g, -g) 



(24) 



-2 . 




FIG. 6: Using the flow equation for S 



(4) 



fa +7i 



FIG. 7: Eq. ggj in its final form. 



r 



Using the relation S ( '{q)A q = 1 (cf. fig. 0, one of 
the two terms on the r.h.s. has been simplified. We can 
see that it just cancels the four-point S contribution in 
eq. 123|) . The remaining one will cancel when eqs. (|21[1 . 
(|22|l are solved for fa. Eq. (|24() is shown in fig.El while 



eq. i|22[l in its final form is displayed as in fig. [7| 

We can diagrammatically process eq. 112 1^ in pretty 
much the same way: integrating by parts and making use 
of the flow equation for the six-point effective coupling, 
we get the diagrams shown in fig. [SJ 



J 



r 



Having a closer look at fig. 03 we see that the terms con- 
taining six-point S vertices cancel out (last diagram on 
the third row and first on the fourth). Moreover, the 
second and third diagrams on the fourth row pair up to 
give a II four-point vertex, which in turn cancels against 
the first diagram on the third row. This is clear once 
this latter is expressed in terms of tree-level vertices. In 
more detail, the one-loop flow equation for the two-point 
coupling at null momentum is just given by fig. except 
with p = and without the C(p 2 ) restriction. With the 
same manipulations we get fig. El again without the or- 
der p 2 restriction and with p = (causing the last term 
in fig. to vanish). Thus we are only left with the total 
derivative in fig. [S] which can be integrated immediately 
to 



s[ 2 \o) 



A„^ 4) (0,0,g, 



(25) 



with no integration constant since in a massless theory 
there must be no other explicit scale apart from the effec- 
tive cutoff. The total derivative term in fig. Q3] vanishes as 
there is no obstruction to exchanging the order of deriva- 
tive and integral signs. 



We still have to deal with the last three diagrams in 
fig- El Of these, we can only "attack" the one with no de- 
pendence upon S, or else we would need to deal with the 
(generic) derivatives of the seed action vertices. So, in- 
tegrating by parts and using the classical flow equations, 
the diagrams shown in fig. 1101 are arrived at. 




FIG. 10: Further processing eq. Hit , from fig- IHI AH external legs carry zero momentum. To integrate by parts the second 
diagram, we use the fact that A q A q = |(Ag)*. 
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FIG. 11: Diagrammatic representation of f3\. 



Collecting our results and solving for f3i we get just the 
five diagrams in fig. 1111 of which only the third remains 
after the renormalization conditions are taken into ac- 



count. In fact, as s£°(0) = S^'ip) 

braic expressions corresponding to the other four dia- 
grams cancel out in pairs. 
In formulae, 

Pi = \ J ^Aa A {c g ^ 4) (0,0, g ,-g)} 2 



(2), 



1, the alge- 



3 n 4 



2 (2^) 4 Jo 
3 



16tt 2 



(26) 



which is the standard one-loop result |32|. 



Moving on to the actual calculation, we start off by 
writing the RG equations for S^iO) and S^\p) , 

p2 

which reduce to purely algebraic relations when the 
renormalization conditions are imposed. They take the 
form: 



fo+ 2 72 - 4Si 2) (0) A 11^ (0) - (0) A S£> (0) 



(4), 



.-"(4) I 



4S{ 2 \0)A Sr'® -§y A q E^>(0, q,-q), (27) 



1 



-.(6)/ 



dqd q |c 9 ^ 4) (0,0,g, -q)\ P2 +72 



-2S { 2 \p)\s¥\p) 



(2), 



(28) 



ED 



and are represented diagrammatically as in figs. 1121 



IV. TWO-LOOP BETA FUNCTION 

This section is devoted to computing the two-loop beta 
function with a generalised seed action. The same proce- 
dure as in the previous section will be followed, namely: 

- introduce integrated kernels in the diagrams containing 
effective action vertices only and integrate by parts so 
as to end up with a total A-derivative contribution plus 
terms where the A-derivative acts on the effective action 
vertices (c/. fig-EJ); 

- make use of the flow equation to process these latter 
further (cf. fig. EI); 

- use the relation between the tree-level two-point vertex 
and the integrated kernel, eq. l(2U|> . i.e. fig. El to simplify 
the structure of the diagrams. At this point, many can- 
cellations should become evident; 

- repeat the above procedure when dealing with the di- 
agrams generated by the use of the flow equations, until 
there is no dependence at all upon the (generic) seed ac- 
tion. 

The one- loop relations, eqs. i|21[l . I|22|) . as well as eq. i|25|) 
will also be used to recast some terms in a more conve- 
nient form. 



02 + 2 72 =A -^2y 





FIG. 12: Graphical representation of eq. 12711 . All the external 
legs have null momentum. 



#2+72 =-2^0j-! ' 

~\ (• ( ^1/ 



2 v? + pJ i 



-X p 



FIG. 13: Graphical representation of eq. 1281 1. 

It is easy to see that the algebraic expressions for the 
diagrams containing the seed action vertices at the high- 
est possible loop order, i.e. the second in fig. 1121 and the 
first in fig. 1131 cancel out when solving for f3% once the 
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renormalization conditions, eq. I|16|) . have been taken into 
account [cf. the comment below fig. I10| . Therefore, in 
what follows, those two diagrams will not appear in the 
graphical representations ofeqs. (|27|) . I|28l) . 

As in the one-loop case, we will start with the easier 
part, i.e. eq. I|28jl. While the first two terms in fig. 1131 
need just expanding at order p 2 , the third has to be fully 
processed. 

Introducing the integrated kernel and making use of 
the flow equation for s[ 4 ^ (p, —p, q, —q) we get the contri- 
butions in fig. El Of these, we can only further simplify 



v and vi, by unfolding the four- and six-point tree-level 
couplings together with the two-point one-loop vertex. 
The way forward is essentially the same as before, but 
with one caveat. Integrating by parts, that is trading 
A q ~ 1/A 2 for A q ~ l/q 2 , might cause the appearance 
of infrared divergences, absent in the first place as the 
ERG equation is well defined in the infrared. Of course, 
these divergences cancel out when all the diagrams on 
the r.h.s. are taken into account, but individual contri- 
butions might be infrared divergent. As an example, let 
us go back to fig. ITU and focus on v [33l| . 



02 + 72 = 




FIG. 14: Unfolding the one-loop four-point effective coupling. 



Upon integration by parts we obtain the contributions 
shown in fig. El 

Of these, the last two are clearly divergent. However 
their divergences must cancel out by construction, as 
the diagram on the l.h.s. of fig. 1151 is infrared finite. In 
the case of the previous example, we could have written 
S[ 2 \q) = [Sf \q)} R + S[ 2} (0), which defines the reduced 
vertex, [Sj 2 ^)]^ The reduced vertex is at least 0(q 4 ), 



and this is what we could have processed further. All the 
diagrams obtained by integration by parts would have 
been finite, containing either a reduced vertex or a A. 
For simplicity's sake, we will continue with the original 
method, bearing in mind that all infrared divergences in- 
deed cancel out, which we will verify at the end of the 
calculation. 



I 

Returning to fig. 1141 and iterating our procedure, we obtain the diagrams in fig. 1161 
I 

I 



Of those, i and iii vanish identically as they are 



(ultraviolet- and) infrared-finite dimensionless integrals, 
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FIG. 15: Appearance of infrared divergences. External legs respectively carry p and —p. 



(3-2 + 72 






01 10 



IV. 



+ 



+W (o 




VI. 
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FIG. 16: Eq. (gHJ in its final form. External legs carry momenta p and — p respectively. The r.h.s. must be evaluated at 0(p 



so do ii and iv after their infrared-divergent parts have 
been cancelled against each other by means of fig. [cf. 
eq. I|25|l]. (N.B. Neither ii nor iv is infrared divergent 
after the A-derivative is taken, however they each yield a 
finite non-zero answer. The strategy is that by combining 
the two diagrams we get one that is finite before differen- 
tiation, and thus zero afterwards.) We are then left with 
v, the wavefunction renormalization contribution to the 
beta function; vi, which will cancel a logarithmic diver- 
gence under a total A9a, and with the last five diagrams, 

(2) 

vii - xi, which all contain Sq (p). Their corresponding 
algebraic expressions will cancel against similar contribu- 
tions from eq. (|27l) by the renormalization conditions, in 
exactly the same fashion as described in fig. ^2 Diagram 



vi is arrived at by combining the expansion of i, iii, iv, 
viii in fig. [21 at order p 2 together with that of other two 
diagrams coming from vi in the same figure and using 

fig- El 

Let us now deal with eq. (|27|) . In complete analogy 
with the one-loop calculation, we first solve the flow equa- 
tion for the two-point coupling at null momentum. The 
two-loop analogue of eq. I|25|) is shown in fig. El (As 
before, the first two diagrams are actually infrared di- 
vergent, but the divergent parts cancel out if eq. I|25|) is 
taken into account.) 

Next, we process the second diagram in fig. ^] After 
several iterations, we arrive at the final form for eq. 127fl . 
shown in fig. 
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FIG. 17: Diagrammatical solution of the flow equation for (0). 




xvi. xvii. xviii. 

FIG. 18: Eq. in its final form. External legs carry zero momenta. 



Again, there are many cancellations. In detail, i, iii, iv mensionless 34], and so do ii and xiii after their infrared- 
and vi vanish identically as they are convergent and di- divergent parts have been cancelled against each other. 
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Moreover, the infrared-divergent part of v results in the 
finite contribution — j3i J, Sq 6 ^ (0, k, — k) which can- 
cels xi exactly. In a somewhat similar fashion, we see 
that the divergent part :35] of viii is equal and opposite 
to that of vii, up to a term that will cancel vi in fig. 1161 
i.e. 

( viii )div = "( vii )div + ft / ^\q, ~Q> k, -k) A fe . 

(29) 

( 2") 

(In more detail, by the renormalization conditions S[ (q) 



has no 0(q 2 ) part, its zeroth order part being that of 
fig. El On the contrary, the 0(q 2 ) part in viii remains 
uncancelled.) Diagrams xi and xii again result from pair- 
ing up several topologically different diagrams and using 
ng.0 



As for the last five diagrams, xiv - xviii, they respec- 
tively cancel the last five terms in fig. ^| when the two 
equations are solved for /?2- Hence, the only diagrams to 
be evaluated are those in fig. HD1 




The first three are already expressed as total deriva- of eq. (|26|l (c/. fig. 120(1 . 
fives, whereas the fourth can be turned into one by means 




FIG. 20: Turning last term in fig. 1191 into a total derivative. 



Although trading a numerical coefficient, 0i, for a di- 
agram might be regarded as an added complication, it 
actually greatly simplifies the rest of the calculation, as 
we now just need to extract the infrared part of each 
diagram. In order to do this, the reduced four-point ver- 
tices are introduced, which by construction (and Lorentz 
invariance) vanish at least as the following powers of mo- 



mentum: 

s R (k) = ^ 4) (o, o, k, -k) - 4 4) (d) = o(k 2 ), 

S R (k, q) = (k, -k, q, -q) - (k, -k, 0, 0) 

-Si 4) (0,0,q,-q)+Sf\d) = O(q-k). 

Eliminating all four-point vertices in favour of the 
reduced ones, we can easily single out the infrared- 
divergent parts, discarding finite contributions which 
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would vanish after the derivative with respect to A is 
taken. 

The third term in fig. 1191 for example, can be written 

as 



^ 4) (0,0,g,-g) A 2 q S { Q 4 \q,-q,k,-k) 

k,q 

xA 2 fe ^ 4) (0,0,fc,-fc) 

= \ j (Sr (q) + l) A 2 (Sr (k, q) + Sr (k) 
+Sr (q) + l) A 2 (Sr (k) + 1 



k.q 



k,q 



2(SR(q) + 
2S 2 (q) - 1 



1 



A 2 q A 2 



A q Afc, 



where S(q) stands for (the un-reduced) 

S^\o,0,q,-q).^ 

The diagram in fig. 1201 has exactly the same structure, 
but with — | in front. The second term in fig. ^3 can be 
simplified to 



3 / \s(q) 2 Al A q A k+q 

Jk,q L 



(30) 



Putting the three together yields 



3 jf ^ [S(q) 2 A 2 A q ( A k+q A q ) ] ' + ~ jf 



k,q 



A 2 . A 2 



k,q 



k . q 



Al A„ ( Atj_„ — A, 



k , (/ 



A 2 . A 2 



A 2 A g (A fc+9 --A 9 



(31) 

where we again used S(q) — Sii(q) + 1 to discard the 
contributions that are infrared finite before differentia- 
tion with respect to A. As regards the first diagram in 
fig. ^3 it amounts to 



1 



k,q 



1 



k.q 



1 

3 Jk,q 



(<So 4) (~p, k,q + p,-k-qfj A fe A q+k A q+p 
(^ 4) (0, k, q, -k - qfj 2 A k A q+k A q+P 
A k A q+k A q+P 



(32) 

as the only non- vanishing contribution comes from taking 
the order p 2 from one of the propagators. 

The integrals H31|M32I) . representing the usual contri- 
butions to the two-loop beta function, can be easily com- 
puted for a generic cutoff function if Bonini's lead [2lJ is 
followed. It consists in "eliminating" one of the cutoff 



functions by writing, say, c q = (c g — 1) + 1 and neglect- 
ing the first term as it is already of order q 2 . This, of 
course, can only be done if the integrand remains UV 
regulated afterwards and if the 0(q 2 ) part indeed makes 
the integral infrared convergent. 

Taking eq. I|32|) as an example, we can trade c k for 1, 
e.g. 



A k A q+k A q+P 



1 

3 Jk.q 



k , q 



l 



Ag_|_fc Ag +p 



1 



Az 



y\2 q+p 



where the last equality follows from taking the derivative 
with respect to A and shifting k — * k — q. 

Averaging over the angles and taking the order p 2 , we 
arrive at the final result (the calculation is detailed in the 
appendix) , 



2 f 1 ( 2_ , 

3j k , q (k-q) 2 Afe U 2Cg+P 



1 / 1 



3 V 16tt 2 



(33) 



The integral in eq. 1)31(1 can be dealt with in pretty 
much the same way (for the details see the appendix): 



As, 9 



A 2 .A 9 (A fe+g --A 9 



= 3 



k.q 



(4c 2 y 

k^q 2 



{k + qf 



1 



= -6 



16tt 2 



and summing up the two contributions the standard re- 
sult follows: 



/?2 
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y 



i 



16ti- 



(34) 
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APPENDIX A: THE EVALUATION OF 
TWO-LOOP INTEGRALS 

This appendix is devoted to the detailed calculation of 

(EH, (El- 
Let us consider first the wavefunction renormalization 
contribution, eq. I|32|l . Rewriting c k as (c k — 1) + 1, we 
see that {c k — 1) already makes the integral convergent in 
the infrared, thus giving a vanishing contribution when 
the derivative with respect to A is taken. Therefore we 
only retain the 1: 



k.q 



A fe A„ +fe A 



q+p 



k . q 



1 



(Al) 
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Taking the A-derivative and bearing in mind the permu- 
tation symmetry in the three momenta: 



(A2) 



where we have also shifted k to k—q. Defining dk to be the 
angle between the Euclidean 4- vectors [36j k and q, i.e. 
k ■ q = kq cos Q k , we can perform the angular integration 
in k to get 



3A2 



1 



4$h 
3A 2 



dkk 3 A k 
Jo J (2tt) 4 (fc - q) 2 



Q+P 



P JO 



dk 



kc k 



max{fc 2 , q 2 } ' 



where fk is the four-dimensional solid angle divided by 
(27r) 4 . Expanding c' q+p to take the order p 2 and defining 
9 q as the angle between q and p, we can integrate over 
the solid angle, d£l q , 



3A 4 



dqq 3 



2A 2 q 



q dk k -^ 



dk 



We now introduce two dimcnsionless variables. 
TP i V = A? and recast the above as 



k 
(A3) 

a; = 



(A4) 



with ci n being c^(a;) and similarly for y. 

Exchanging the order of integration, so that the inte- 
gral over y is performed first, and temporarily discarding 
the numerical factor in front of the integral 



dx < 
1 

1 



dy 



V_ r w 
2 y 



z ./o Jo 



, c x X , 



Hence, eq. (13211 amounts to 



Iff 
4 3 



1 



1 



3 V 16tt 2 



(A5) 



As far as eq. (|31J1 is concerned, we can follow the same 
strategy and write c k +q = {ck+q —c q ) + c q , retaining only 
c„: 



= 3 



A?.A 9 ( 



1 . 

^k+q ~ 2 9 



1 



k 4 q 2 y(k + q) 2 



1 

V 



Averaging over the angles and rewriting in terms of 
x, y we get 




(c 2 x y 



dx / dy 



dx / dy 



(A6) 



Again exchanging the order of integration and using 
the fact that the integrand is invariant under the ex- 
change x «-> y, we see the last two lines in eq. l|A6ll are 
equal and opposite, while the first can be rewritten as 



c 

X Jo 



- I dyc 2 y 



3 r 2 
= - hm ifiS. 



dyc 2 y 



( 1 



V16tt 2 
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